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Abstract 

We show that the assumption of quasiperiodic boundary conditions (those that 
interpolate continuously periodic and antiperiodic conditions) in order to compute 
partition functions of relativistic particles in 2+1 space-time can be related with 
anyonic physics. In particular, in the low temperature limit, our result leads to the 
well known second virial coefficient for anyons. Besides, we also obtain the high 
temperature limit as well as the full temperature dependence of this coefficient. 
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In nature one encounters particles which obey Bose-Einstein (BE) or Fermi-Dirac 
(FD) statistics. However, speculations concerning the existence of particles which obey 
intermediate (fractional) statistics between BE and FD have been raised for a long time. 
In recent years this kind of conjecture has gained force with its possible relation with 
condensed matter effects like high-Tc superconductivity and the fractional quantum Hall 
effect §. 

In (2 + 1) dimensional field theory, fractional statistics can be obtained by including 
a Chern-Simons term |J in the Lagrangian density describing the interaction of matter 
with gauge fields. Particles with intermediate statistics between the BE and FD cases 
in two space dimensions are usually called anyons. Their properties can be understood 
from the braid group or equivalently assuming that they are hard core indistinguishable 
particles in two space dimensions so that the configuration space becomes a manifold Ai 
with a non-trivial topology in the sense that Ui(Ai) = for N particles, where 
denotes the braid group with N generators. Hence, it can be shown that when anyons 
turn around each other the corresponding wave function (or the Feynman propagator for 
the system) acquires a non-trivial geometric phase An inherent difficulty of the 

anyon field description is its non-locality 0. This makes it difficult to calculate beyond 
the second virial coefficient || even in the case of "free" anyons, contrary to what happens 
for free bosonic and fermionic gases for which the full partition function is known. 

Here, we shall discuss a different approach to the interpolation between bosons and 
fermions. In this paper we shall compute a generalized partition function which contains 
the relativistic partition functions for bosons and fermions as particular cases, and also 
interpolates continuously between these cases. From this partition function we obtain 
the second virial coefficient and show that, in the low energy (temperature) limit it re- 
produces correctly the results known for anyons in the literature ||. We also find the 
high temperature and the full temperature dependence of the relativistic second virial 
coefficient. 

Let us briefly review how to obtain the partition function for ideal relativistic bosonic 
and fermionic gases, from which we will generalize to the anyonic case. The partition 
function for a system described by a Hamiltonian H with chemical potential \x can be 
written as 

Z = exp{-(3tt} = Tr e -^ +A ^), (1) 

where Q is the free energy and TV is the conserved charge of the system. As it is well known 
for relativistic charged massive bosonic fields one can express this partition function as a 
determinant PlflT0|, namely, 

Z = [det(-£> 2 + M 2 ) | ] ~* ; (bosons) , (2) 

where D 2 is the square of the covariant derivative, D u , including the chemical potential, 
D u = (do+i/j,, di) and M is the mass of the field. This prescription for the inclusion of the 
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chemical potential \x coincides with the one given by eq. (0]) except for the introduction 
of an imaginary part to Z, which is not physically relevant. So unless otherwise specified, 
we are considering that the real part of Z has been taken. Note also that, the chemical 
potential appears squared in the relevant operator — D 2 + M 2 . This is a consequence of 
representing the partition function in phase space and then integrating over momenta (see 



1Q| for details). The label P means that the eigenvalues of this operator are subjected to 



periodic boundary conditions and hence are given by 

Kk = K + i^f + k 2 + M 2 , (3) 



where ui n = 2rnr/{3, with n G 2, are the Matsubara frequencies [|TT| for bosonic fields and 
k e R N . 

As for the case of bosons, one can also write the partition function for free relativistic 
fermions as a determinant 



Z = det D (0) - M)\ A 
det(-D 2 + M 2 " 



(fermions), (4) 



where the determinants are calculated over the corresponding operator and detp means 
the inclusion of the calculation over Dirac indices (d is the dimension of the Dirac rep- 
resentation which hereafter we take d = 2). The subscript A means that the eigenvalues 
of — D 2 + M 2 are now computed with antiperiodic boundary conditions. Hence, they are 
given by (H), but now the Matsubara frequencies are ui n = (2ir/ (3){n + 1/2), due to the 
antiperiodic boundary condition. 

The remarkable fact about the (last line of) equation (|j) is the well known but few 
explored character that it states that it is possible to compute the partition function for 
fermions through an equation where spin does not appear at all. The information that 
the partition function (f|) corresponds to fermions is solely described by the fact that 
the determinant of the Klein-Gordon operator (-D 2 + M 2 ) is computed over antiperiodic 
eigenf unctions. Naturally, this antiperiodicity comes from the fact that fermions anticom- 
mute when they are exchanged in space, i. e., the wave function of the system acquires a 
sign (—1). So, on thermodynamical grounds the only necessary information to compute 
a partition function for free fermions is that they obey energy conservation (expressed 
through the Klein-Gordon operator) and that their wave functions are antiperiodic in the 
Euclidean time. 

Now, we want to extend this fundamental result to the case of (quasi) particles that 
obey an intermediate statistics between the bosonic and fermionic cases. As is well known 
these particles live in 2+1 dimensional space-time and they can be described by specific 
equations of motion like the one of Jackiw and Nair [|l2| . This equation of motion governs 
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the dynamics of that particles and contains detailed information on their algebraic prop- 
erties, including the braid group as well. However, the task of computing a determinant 
of an operator obtained from such equation is very involved and we are not going in this 
direction in this work. 

An alternate route, that we are going to follow here, is to use very few but essential 
information which are necessary for a thermodynamical description of a system of such 
particles. These informations are only three: first these (relativistic) particles should 
obey energy conservation and so in spite that they are described by some involved specific 
2+1 dimensional equation they should also obey the Klein-Gordon equation; second the 
determinant of the Klein-Gordon operator should be calculated using some appropriate 
boundary condition; Finally the partition function is defined by a certain power a (to 
be discussed later) of the determinant of the Klein-Gordon operator, under a chosen 
boundary condition. Of course, if the boundary condition is periodic in Euclidean time 
and a = — 1 we will be describing bosons as is expressed in equation (Q). On the other 
hand, if we impose antiperiodic boundary conditions and a = +1 we will be describing 
fermions as in equation (|j). As we want to describe particles that obey some statistics 
that continuously interpolate the bosonic and fermionic cases, it seems rather natural 
to choose a boundary condition that interpolates these two well known cases. Such a 
condition, which we call quasiper iodic, can be written as 



where 9 is a continuous parameter defined in the interval [0, tt}. Note that putting 9 = 
we reduce the above condition to a periodic boundary condition and if 9 = it it becomes 
an antiperiodic one. 

The power a of the determinant is introduced to fit the bosonic and fermionic cases. 
As is well known in these cases the powers —1 and +1, respectively, come from the 
definition of the measure over c- number (bosons) and Grasman variables (fermions). For 
the anyonic case it would be necessary to define a set of variables which could interpolate 
the properties of c-numbers and Grasman variables. This approach seems to be related to 
q-deformed calculations, but here we are just going to leave the power a of the determinant 
free and check, at the end of the calculations that both o = ±1 will be related to anyonic 
physics. This may be explained remembering that conventional anyons are constructed 
from fermions or bosons (with a fixed built in measure) adding a Chern-Simons term. 
Later on, we will also show that keeping 9 fixed and varying a we could describe as well 
the interpolation between bosons and fermions. 

Note that the quasiperiodic condition (||) can also follow from an analogy of the spatial 
behavior of bosons, fermions and anyons. This behavior is the well known property that 
bosonic wave functions are symmetrical while fermionic ones are antisymmetrical and 
anyonic wave functions acquire a phase exp(i9) when these particles are turned around 




i9 



(5) 
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each other [Q, [13|]. Going to the finite temperature case, these properties imply that 
bosonic fields are periodic in Euclidean time while fermionic are antiperiodic. This analogy 
is not complete for anyons and we are assuming here that an intermediate continuous 
condition can explain at least some of the thermodynamical properties of anyons, as we 
are going to show bellow. 

Then, we are not proving that anyons obey the quasiperiodic boundary condition, eq. 
(|5]), and we are not able at this point to do so, but we are using a very appealing analogy 
to support it. The condition @ is our basic assumption from which we are going to show 
that the partition function that we will calculate can be mapped onto the corresponding 
results for conventional anyons. In fact we do this for the only exact result known for 
anyons in this matter which is the second virial coefficient obtained by Arovas et. al. ||. 
The calculation of higher virial coefficients and its comparison with the perturbative or 
numerical ones known in the literature is done in a forthcoming work [|TJ] . 



The price we have to pay in our approach is that we are loosing detailed quantum 
mechanical information as those associated with the braid group which become very 
involved when many particles come into play. The obvious advantage of this global 
approach is that it keeps only the essential thermodynamical information which suffices 
for calculating the virial coefficients, for example. 

Note also that, in our approach there is no reason for restricting the space-time dimen- 
sions, despite that anyons live only in two space dimensions. This is not surprising since 
other approaches towards generalized statistics, like the Haldane's generalized exclusion 
principle []15|| , can be formulated in spaces with arbitrary number of dimensions. In fact, 
at the end of our calculations we shall restrict them to two space dimensions in order to 
show the equivalence of our approach with the conventional one for anyons. 

Let us now compute the following determinant of the Klein-Gordon operator: 



Z 



det(-D 2 + M 2 



(6) 



where the subscript 9 means that the quasiperiodic condition (||) is assumed and we 
introduced the parameter a to be able to reproduce correctly the bosonic and fermionic 
particular cases. Observe that when 9 = and a = — 1 we reobtain the bosonic partition 
function, while for 9 = tt and o = +1 we have the fermionic case. For generality, we are 
going to calculate the above determinant in N + 1 space-time dimensions. The eigenvalues 
of the Klein-Gordon operator in this case are: 



A 



nkd 



2mr 9 



+ k 2 + M 2 . (7) 



The determinant (^|) is a generalization of its quantum mechanical (0 + 1) dimensional 
version (with fi = 0) which has been calculated using Green functions |16| and the zeta 
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function method [17]. Here, we are going to calculate this determinant also using the 
generalized zeta function method through the basic formula [ITH 



det A = exp < — — £(s; A) 



s=0 



where the generalized zeta function is defined by £(s; A) = Tr A~ s and an analytical 
continuation of £(s; A) to the whole complex plane of s is tacitly assumed. 

For the case at hand, with A = —D 2 + M 2 under ^-periodic boundary conditions, the 
eingenvalues being given by eq. fl7|), the generalized zeta function reads 



Va 



N 



(2*0 



N 



E 



dk k 



N-l 



2mr 6 



+ k 2 + M 2 



(9) 



where V is the volume, a at is the area of the unit hypersphere, both in N space dimensions. 
The integral is expressible in terms of the Beta function, and we obtain 



2 r(*) V + 



where we have defined 



C 



N 



Va N 

N ' 



(2tt 



/3M 
~2tV 



ipH + 9 
2tt 



(10) 



The sum appearing in eq. (fTD[) , well defined for 9ft(s) > 1 + —, is a generalization of the 
usual Epstein function KS| but as one can see from formula (^) , that the region of interest 
contains the point s = 0. The analytic continuation of this sum for the whole complex 



s-plane can be written as [20 



E W 2 + (n + X ) 2 ]-^ 
'T(s 



71 



r(s- 



N ■ 



N 
2 



+oo 



2s-N-l 



Aj2cos(27inx).{ 



— — — / ^ 

2 2 K s _N__i{2irnv) 



n=l 



, (12) 



where K a (z) is the modified Bessel function of order a. Taking the derivative of ((s,A) 
with respect to s, using that lim s ^ Q [l/T(s)] = and lim s _ [(d/ds)T(s)] = 1 we have, 
apart from an irrelevant term which is linear in /3 and independent of [i: 



s=0 



2V^C N T(^) N g cos(2vrn X ) (™ 



-i(iV+l) 



2 



(JV+1) (27rnv) 



(13) 
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So, using the fact that lndet^l = —(d/ds)((s,A)\ s= o, we can find the generalized deter- 
minant (|G|) which interpolates continuously the partition functions, or equivalently the 
free energy, eq. ([I]), for relativistic boson and fermion gases with chemical potential /i in 
(N + l)-dimensions: 

n((3,ii) = -4 In 2 = ~<r In det A 

P P 

Va N 1 / , 2M\* (W+1 > .N. 

+oo 1 

x ^cos(n^)cosh(n/?/i)(-)5( 7V+1 )iri (7V+1) (n/?M) , (14) 

n=l n 2 

where Cjy, v and x have been taken from eq. (|TT|). Note that in the above formula, 
only the real part of the free energy f2(/3, /i) was taken into account, according to the 
prescription of introducing the chemical potential as an imaginary time-component gauge 



potential [10|. If we particularize the parameters a and 9 to the bosonic (a = — 1 and 
9 = 0) and fermionic (er = +1 and 9 = 7r) cases we shall find precisely the results known 
in the literature pTl-pj. 



To show the equivalence of this partition function with the conventional theory of 
nonrelativistic anyons, we first recall the cluster expansion 

^hi£ = £ k(V,T)z l , (15) 

where z = exp is the fugacity and bi(V, T) are the usually called cluster coefficients. As 
Z = exp(— (3Q), we substitute the expression of /x), eq. (p!4|), into the above equation 
and find that the cluster coefficients, in our case, are given by 

JV+l 

/ M \ ~2- x_ N 

b±JV,T) = - 2a cos(n0 )KN+i(n(3M) ; (n = 1, 2, 3, ...) . (16) 

V 27m / 2 

The ± sign for the cluster coefficients are related to particles and antiparticles. Note 
that, in a conventional nonrelativistic system only particles (or antiparticles) are present, 
so in that case we would have found only one set of natural cluster coefficients. As they 
are in fact equal, as required by particle-antiparticle symmetry, this will not modify our 
analysis. 

Now, using the standard formula which relates the second virial coefficient with the 



cluster ones [23| 
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we have in N space dimensions the exact expression for the relativistic second virial 
coefficient interpolating the bosonic and fermionic cases: 



B ^ = To Km) 



N+l 

^ 2 (3^(1 



tan 



Kn+i (2/3M) 



Kn+i(pm) 



To write the ratio of the two Bessel functions in a more familiar form we can use the 
asymptotic expression for them 



K v [x) = \l — e 
2x 



Av 2 -l (4z/ 2 - l)(4z/ 2 -3 2 ) 



8:r 



2! (8x) 



(19) 



valid when — n/2 < aigx < tt/2. Once we are interested in relating our relativistic 
partition function to the conventional formulation of nonrelativistic anyons, we must 
reduce our formula by taking the low temperature (nonrelativistic) limit (3M » 1 on 
equation ( |i~8|) of the second virial coefficient. In this case, the above asymptotic expression 
for the Bessel functions can be greatly simplified to 



K »( x ) - V y x e x ; 



(x » 1) 



(20) 



Form now on, we will restrict our analysis to the particular case of interest, which is 
N = 2 space dimensions where anyons live. So considering the low temperature approxi- 
mation and using eq. (p0|), we find that eq. ( p~8|) reduces to (h = c = 1): 

1 7T(3 



B(T) 



2a M 



[l-t&n 2 9) 



((3M » 1). 



(21) 



We show now that this is equivalent to the results of Arovas et. al. for the second virial 
coefficient. If one starts from charged fermions interacting with a magnetic flux tube with 
intensity (f> = ahc/e, one gets [H] 



B{T) = \\l (1-25 2 ) 

or analogously, starting from bosons the result is 

1 



B{T) = --X 2 T (1 - A\5\ 



25 2 



(22) 



(23) 



where A 2 - = l?^ is the thermal wavelength and 5 is the non-integer part of the statistical 
parameter a, i. e., a = 2j + 1 + 5 for fermions, eq. (0), and a = 2j + 5 for bosons, eq. 
3D, j = 0, ±1, ±2, ... in both cases. 
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To compare these results, we note that, in our case, if we start from fermions we 
should put cr = +1, identify tan 2 9 = 25 2 and get the same answer for the second virial 
coefficient of anyons, as in eq. (|22|). Note that these fermions came about through the 
choice 9 = ii, but we can redefine 9 = 9' + nn such that — n/2 < 9' < +n/2, since the 
9 dependence appears only in the second virial coefficient through tan 2 9. In fact, if we 
want to interpolate the virial coefficients of bosons (— A^/4) and fermions (+Af./4) we 
must restrict tan 2 9 < 2. This restriction, which can be written as — arctan \[2 < 9' < 
arctan a/2, is exactly the one which is needed to map our second virial coefficient onto the 
one of Arovas et. al. 

In the other case, if we start with bosons, we put o = — 1 and identify tan 2 # = 
2|5|(2 — \5\), in a complete analogous situation where the same restrictions on 9, from 
the fermionic case, applies here as well. One should also note that, the Arovas et. al. 
second virial coefficient presents cusp singularities for Act = 2, which are not present in 
our result since the coefficient we found is analytical in its parameters except for poles 
at 9 = (n + l/2)7r, where n = 0, ±1, ±2, ... This difference comes from the fact that our 
result also contains values for B(T) which extrapolate the interval (— +A^/4). 

A simpler situation comes up if we consider the particular cases of small angles 9 ~ 
or 9 ~ 7r, so that 

tan 2 9 ~ 9 2 . 

In these situations the relation between the quasiperiodic 9 and the statistical parameter 
a (or 5) becomes simply (considering the fermionic anyon model given by (|22|)): 

V2a = 9. 

So, in the region where anyons are close to fermions (or bosons), where usually per- 
turbation calculations are done for other quantities like the higher virial coefficients, 9 
plays the role of the statistical parameter a up to a constant factor \/2. Of course, the 
exact relation between these parameters is nonlinear as found above, but in that case the 
analysis is much more involved and we do not have until now a deeper understanding of 
this nonlinear behavior. 

To understand how we succeed in describing the statistical mechanics of anyons in 
this simple way, we may say that, the interpolating parameter 9 discussed here plays the 
role of a topological constant gauge field = (A , 0), since the ^-boundary condition for 
the Klein-Gordon operator, both in the bosonic as well as in the fermionic case, implied 
that the Matsubara frequencies change according to u n — > uj n + 9/ j3, which can be viewed 
as a shift in the time derivative operator, <% — > do + i9/(3. As the time coordinate, 
xq, is compactified to the interval (0,f3), this introduces the non-trivial topology usually 
associated with anyons but in a different and unexpected manner. This prescription, which 
works only for the finite temperature case, allows us to relate directly the quasiperiodic 
boundary condition with intermediate statistics, or anyon description. 



9 



Another interesting point in our formulation is that we can interpolate the bosonic and 
fermionic second virial coefficients, from two different points of view. First, keeping a fixed 
to ±1 and varying 9 we were able to describe this interpolation very close to the anyon 
description starting from bosons and fermions with their original equations of motion, or 
equivalently the bosonic and fermionic determinants, with the inclusion of a topological 
term, in our case the ^-boundary condition, as we described above. Alternatively, we can 
keep 9 fixed to or 7r and vary continuously a between +1 and —1, so that we should 
identify cr -1 = 1 — 2S 2 in the fermionic case and a -1 = 1 — 2|5|(2 — \S\) in the bosonic 
case. This seems to be more close to a q-deformed calculation but this connection is, at 
least for the moment, far from being trivial. 

We can go even further in the determination of the relativistic anyon second virial 
coefficient, since we do not need to use the low temperature approximation in our calcu- 
lations. This is so because, in our exact result, eq. flTB|), we can also use (in N = 2) the 
exact expression for the Bessel function of order 3/2: 

Ks(x) = ,[Ze~ x (l + -) . (24) 



2 V 2x V x 

Substituting this expression in the second virial coefficient, eq. (|IBD , instead of the 
approximate one, eq. (|20|) , we find the exact relativistic result for the anyon second virial 
coefficient: ^ 

B(T) = — ^(1- tan 2 9) + + — ) . (25) 

y J 2a M ' \ 2(3M ) \ (3M ) y ! 

Note that the possibility of finding this exact result for the second virial coefficient is 
a peculiarity of odd-dimensional space-times, since for these cases the relevant Bessel 
functions are of half-integer order for which the asymptotic expansion became exact (see 
eq. flnp). Naturally, the low temperature result, eq. fl2"T|), can be easily obtained from 
fl25|) just taking the low temperature limit (3M » 1. 



Furthermore, we can also find the high temperature limit (extremely relativistic case) 
of expression ( p5|) taking (3M « 1 so that we find 

BIT) = — 7T/? 2 (1 - tan 2 9) ; (BM « 1) . (26) 

Note that the statistical correction (1 — tan 2 9) is independent of the high/low tem- 
perature regime, but the mass dependence disappears in the high temperature case, while 
the power of temperature just changes from (3 to /3 2 . 

As a final comment let us mention that it seems interesting to investigate the extension 
of these results to other related situations. One of these is the calculation of the higher 
virial coefficients which are very promising since they can be exactly evaluated in our 
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approach |14| . Another situation of interest is the inclusion of an external magnetic field 
since in most applications of anyonic physics, like the fractional quantum Hall effect, 
strong external fields are present. We should report on this elsewhere. 
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